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ABSTRACT 

The paper experimentally and theoretically studies 
the effects of periodic unsteady wake flow and 
aerodynamic characteristics on boundary layer 
development, separation and re-attachment along the 
suction surface of a low pressure turbine blade. The 
experiments were carried out at Reynolds number of 
110,000 (based on suction surface length and exit 
velocity). For one steady and two different unsteady 
inlet flow conditions with the corresponding passing 
frequencies, interm ittency behavior were experimentally 
and theoretically investigated. The current investigation 
attempts to extend the intermittency unsteady boundary 
layer transition model developed in previously to the 
LPT cases, where separation occurs on the suction 
surface at a low Reynolds number. 

The results of the unsteady boundary layer 
measurements and the intermittency analysis were 
presented in the ensemble -averaged, and contour plot 
forms. The analysis of the boundary layer experimental 
data with the flow separation, confirms the universal 
character of the relative intermittency function which is 
described by a Gausssian function. 
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INTRODUCTION 


In recent years, gas turbine engine aerodynamicists 
have focused their attention on improving the efficiency 
and performance of the low pressure turbine (LPT) 
component. Previous research has shown that a 
reduction in the blade numbers can be achieved without 
substantially sacrificing the efficiency of the LPT 
blading. For aero-engines this reduction contributes to 
an increase in thrust/weight ratio, thus reducing the fuel 
consumption. Contrary to the high pressure turbine 
(HPT) stage group that operates in a relatively high 
Reynolds number environment, based on the operating 
conditions, the LPT experiences an adverse change in 
Reynolds number ranging from 50,000 to 250,000. 
Since the major portion of the boundary layer, 
particularly along the suction surface is laminar, the low 
Reynolds number in conjunction with the local adverse 
pressure gradient makes it susceptible to flow 
separation, thus increasing the complexity of the LPT 
boundary layer aerodynamics. The periodic unsteady 
nature of the incoming flow associated with wakes that 
originate from upstream blades substantially influences 
the boundary layer development including the onset of 
the laminar separation, the extent of the separation 
bubble, and its turbulent re- attachment. Of particular 
relevance in the context of LPT aerodynamics is the 
interaction of the wake flow with the suction surface 
separation bubble. While the phenomenon of the 
unsteady boundary layer development and transition in 
the absence of the separation bubbles has been the 
subject of intensive research that has led to better 
understanding of the transition phenomenon, grasping 
the multiple effects of mutually interacting parameters 
on the LPT boundary layer separation and their physics 
still requires more research. To fully understand the 
basics involving the separation bubble phenomenon, an 
intermittency analysis has been employed to extend the 
intermittency unsteady boundary layer transition model 
developed in [1], [2], [3] to the LPT cases, where 
separation occurs on the suction surface at a low 


Reynolds number. 

Studies by Abu-Ghannam and Shaw [4], Gostelow 
and Blunden [5], and Dullenkopf and Mayle [6], were 
conducted to determine the effect of ffee-stream 
turbulence and pressure gradient on the spot production 
rate and the intermittency factor. Significant 
contributions to steady and unsteady boundary layer 
research was made by Pfeil and his co-researchers ([7], 
[8], [9], [10], [11]). Pfeil and Herbst [12], utilizing the 
squirrel cage-type wake generator and a flat plate, 
developed a wake-induced transition model that is now 
generally accepted as accurate. They also showed that 
the boundary layer grew naturally in between the 
induced transition regions by wakes. Comprehensive 
investigations on the effect of periodic unsteady flow on 
a curved plate were performed by Schobeiri and Radke 
[13], and Schobeiri et al. [14]. They showed that an 
increase in wake passing frequency as a result of 
reducing the wake spacing results in changing the wake 
turbulence structure, and also a shift of transition region 
towards the leading edge. Experiments for the effect of 
unsteady wake flow on the boundary layer transition 
were also conducted by Walker [15], Paxson and Mayle 

[16] , and Orth [8]. Paxson and Mayle investigated the 
effect of unsteady passing wakes on the laminar 
boundary layer near the stagnation region. Dullenkopf 
and Mayle [3] proposed a time averaged transition 
model. Few of these researchers have addressed the 
effect of wake frequency and the structure on the 
boundary layer transition. 

The transition process was investigated by Emmons 

[17] through the turbulent spot production theory. This 
theory was later promoted by Dhawan and Narasimha 

[18] , who found the intermittency factor for natural 
transition. Unlike the steady boundary layer transition 
case, the calculation of intermittency function under the 
unsteady flow situation is a difficult task because of the 
ffee-stream turbulence distribution, which is periodically 
changing from almost non-turbulent to high turbulent 
intensity values. The process ofturbulent/non-turbulent 
decisions from the instantaneous signals measured under 
these unsteady conditions is reviewed by Hedley and 
Keffer [19]. They proposed derivatives of velocity 
signals as the detector function to identify the turbulent 
and non-turbulent parts in the signals. This method was 
also used by Antonia and Bradshaw [20], Kovasznay, et 
al. [21], and Bradshaw and Murlis [22]. Mayle and 
Paxson [13] and Mayle [23] used a similar method for 
unsteady flows. 

Developing an unsteady transition model is essential 
to accurately predict the unsteady boundary layer 
characteristics such as skin friction and heat transfer 
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coefficients. With an appropriate transition model, it is 
possible to numerically solve the boundary layer 
equations using different methods such as those 
proposed by Launder and Spalding [24], Crawford and 
Kays [25], and Schmidt and Patankar [26]. 
Implementing such a model in an unsteady Navier- 
Stokes code enables reliably predicting the 
turbomachinery profile loss coefficients and thus, the 
efficiency. 

Based on the fundamental investigations of the 
velocity and the turbulence structure of the impinging 
wakes and their interaction with the boundary layer, 
Chakka and Schobeiri [1] developed an intermittency 
based unsteady boundary layer transition model. The 
analysis revealed a universal pattern for the relative 
intermittency function for all the frequencies and 
pressure gradients investigated. However, the above 
investigations were not sufficient to draw any 
conclusion with regard to an eventual universal 
character of the relative intermittency function. Further 
detailed investigations of the unsteady boundary layer 
on a high Reynolds number turbine cascade by Schobeiri 
et al. [27], [28] and its subsequent analysis [2] and [3] 
verified the universal character of the relative 
intermittency function. 

The current investigation attempts to extend the 
intermittency unsteady boundary layer transition model 
developed by Schobeiri and his coworkers ([1], [2], [3]) 
to the LPT cases, where a massive separation occurs on 
the suction surface at a low Reynolds number at the 
design and off-design points. Furthermore, the 
experimental results are intended to serve as benchmark 
data for a comparison with numerical computation using 
DNS or RANS-codes. 

EXPERIMENTAL RESEARCH FACILITY 

T o investigate the effect of unsteady wake flow on 
turbine and compressor cascade aerodynamics, 
particularly on unsteady boundary layer transition, a 
multi-purpose large-scale, subsonic research facility was 
designed and has been taken into operation since 1993. 
Since the facility in its original configuration is described 
in [27], [28] and [29] only a brief description of the 
modifications and the main components is given below. 
The research facility consists of a large centrifugal 
compressor, a diffuser, a settling chamber, a nozzle, an 
unsteady wake generator, and a turbine cascade test 
section as shown in Figure 1 . The compressor with a 
volumetric flow rate of 15 m 3 /s is capable of generating 
a maximum mean velocity of 1 00 m/s at the test section 
inlet. The settling chamber consists of five screens and 


one honeycomb flow straightener to control the 
uniformity of the flow. 

Two-dimensional periodic unsteady inlet flow is 
simulated by the translational motion of a wake 
generator (see Figure 2), with a series of cylindrical rods 
attached to two parallel operating timing belts driven by 
an electric motor. To simulate the wake width and 
spacing that stem from the trailing edge of rotor blades, 
the diameter and number of rods can be varied. The rod 
diameter, its distance from the LPT blade leading edge, 
the wake width and the corresponding drag coefficient 
is chosen according to the criteria outlined by Schobeiri 
et al. [30]. The belt-pulley system is driven by an electric 
motor and a frequency controller. The wake-passing 
frequency is monitored by a fiber-optic sensor. The 
sensor also serves as the triggering mechanism for data 
transfer and its initialization, which is required for 
ensemble-averaging. This type of wake generator 
produces clean 2-dimensional wakes, whose turbulence 
structure, decay and development is, to a great extent, 
predictable [30]. The unsteady boundary layer transition 
and heat transfer investigations [1], [2], [27], [28] 
performed on this facility serve as the bench mark data 
for validation of turbulence models, transition models, 
and general code assessments. 

T o account for a high flow deflection of the LPT 
cascade, the entire wake generator and test section unit 
including the traversing system, was modified to allow 
a precise angle adjustment of the cascade relative to the 
incoming flow. This is done by a hydraulic platform, 
which simultaneously lifts and rotates the wake 
generator and test section unit. The unit is then attached 
to the tunnel exit nozzle with an angular accuracy less 
than 0.05°, which is measured electronically. 

The special design of the facility and the length of 
the belts (L^ = 5,000 mm) enables a considerable 
reduction of the measurement time. For the present 
investigation, two clusters of rods with constant 
diameter of 2 mm are attached to the belts as shown in 
Figure 2. The two clusters with spacings S R = 160 mm 
and S R = 80 are separated by a distance which does not 
have any rods, thus simulating steady state case (S R = 
°°). Thus, it is possible to measure sequentially the effect 
of three different spacings at a single boundary layer 
point. T o clearly define the influence domain of each 
individual cluster with the other one, the clusters are 
arranged with a certain distance between each other. 
Using the triggering system mentioned above and a 
continuous data acquisition, the buffer zones between 
the data clusters are clearly visible. 
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® Static pressure blade @ Timing bells, rod attachments ( 9 ) Inlet nozzle 

(2) Blade with hot film sensors Transition duct © Hydraulic cylinders 

© Vake generating rods @ Straight duct © Pivot point 

® Wake generator ® Traversing system © Wake generator e-motor 


Large silence chamber with 
honeycom and five screens 
Telescope supprt 

Honeycomb flow straightener 
Traversing slots 


Figure 1. Turbine cascade research facility with the components and the adjustable 


test section 



Figure 2. Wake Generator 


4 


Copyright © 2005 by ASME 








Table Is Parameters of turbine cascade test section 


Parameters 

Values 

Parameters 

Values 

Inlet velocity 

V fa = 4 m/s 

Inlet turbulence intensity 

Tu in = 1.9 % 

Rod translational speed 

U= 5.0 m/s 

Blade Re-number 

Re =110,000 

Nozzle width 

W = 200.0 mm 

Blade height 

h B = 200 mm 

Blade chord 

c = 203.44 mm 

Cascade solidity 

G= 1.248 

Blade axial chord 

c^ = 182.85mm 

Zweifel coefficient 

t A = 1-254 

Blade suction surface length 

L ss = 270.32 mm 

Cascade angle 

(p = 55° 

Cascade flow coefficient 

$ = 0.80 

Cascade spacing 


Inlet air angle to the cascade 

a, = 0° 

Exit air angle from the cascade 

a 2 = 90° 

Rod diameter 

D r = 2.0 mm 

Rod distance to lead, edge 

L r = 122 mm 

Cluster 1 (no rod, steady) 

S R = 00 mm 

Q - parameter steady case 

Q = 0.0 



Q - parameter for cluster 1 

Q = 1.59 



Q - parameter for cluster 2 

Q = 3.18 


The data analysis program cuts the buffer zones and 
evaluates the data pertaining to each cluster. 
Comprehensi ve preliminary measurements were carried 
out to make sure that the data were exactly identical to 
those, when the entire belt length was attached with 
rods of constant spacing, which corresponded to each 
individual cluster spacing. The cascade test section 
shown in Figure 1, located downstream of the wake 
generator, includes 5 LPT blades with a height of 200.0 
mm and the chord of 203.44 mm. For boundary layer 
investigations, five identical “PakB” airfoils designed by 
Pratt & Whitney were implemented whose cascade 
geometry is given in T able 1 . 

The blade geometry resembles the essential feature 
such as the laminar boundary layer separation that is 
inherent to typical LPT blades. The blade geometry was 
made available to NASA researchers and academia to 
study the specific problems of LPT flow separation, its 
passive and active control and its prevention. As shown 
in [27], this blade number is necessary and sufficient to 
secure a spatial periodicity for the cascade flow. The 
periodicity is demonstrated by the pressure distributions 
of the blade number 2 and 4, shown in Figure 1 These 
blades were specially manufactured for measurement of 
pressure and showed identical pressure distributions. A 


computer controlled traversing system is used to 
measure the inlet velocities and turbulence intensities, as 
well as the boundary layers on suction and pressure 
surfaces. The traversing system is vertically mounted on 
the plexiglass side wall. It consists of a slider and a lead 
screw that is connected to a DC-stepper motor with an 
encoder and decoder. The optical encoder provides a 
continuous feedback to the stepper motor for accurate 
positioning of the probes. The system is capable of 
traversing in small steps up to 2.5 |lm, which is 
specifically required for boundary layer investigations 
where the measurement of the laminar sublayer is of 
particular interest. 

INSTRUMENTATION, DATA ACQUISITION, 
AND DATA REDUCTION 

The data acquisition system is controlled by a 
personal computer that includes a 16 channel, 12-bit 
analog-digital (A/D) board. Time dependent velocity 
signals are obtained by using a commercial 3 -channel 
(TSI), constant temperature hot-wire anemometer 
system that has a signal conditioner with a variable low 
pass filter and adjustable gain. A Prandtl probe, placed 
upstream of the diffuser, monitors the reference velocity 
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at a fixed location. The pneumatic probes are connected 
to high precision differential pressure transducers for 
digital readout. Several calibrated thermocouples are 
placed downstream of the test section to constantly 
monitor the flow temperature. The wake generator 
speed and the passing frequency signals of the rods are 
transmitted by a fiber-optic trigger sensor. The passage 
signals of the rods are detected by the sensor using a 
silver-coated reflective paint on one of the belts. This 
sensor gives an accurate readout of the speed of the 
wake generator and the passing frequency of the rods. 
The signals of the pressure transducers, thermocouples, 
and trigger sensors are transmitted to the A/D board and 
are sampled by the computer. T o ensure the cascade 
periodicity, the second and fourth blades are 
instrumented each with 48 static pressure taps. Two 
adjacent blades are used for boundary layer 
measurement. The taps are connected to a scanivalve, 
which sequentially transferred the pressure signals to 
one of the transducers that was connected to the A/D 
board. 

The unsteady data are taken by calibrated, custom 
designed miniature, single hot wire probes. At each 
boundary layer position, samples were taken at a rate of 
20kHz for each of 100 revolutions of the wake 
generator. The data were ensemble -averaged with 
respect to the rotational period of the wake generator. 
Before final data were taken, the number of samples per 
revolution and the total number of revolutions were 
varied to determine the optimum settings for 
convergence of the ensemble-average. 

For the steady state case, the instantaneous velocity 
components are calculated from the temperature 
compensated instantaneous voltages by using the 
calibration coefficients. The instantaneous velocity can 
be represented in the following form: 

V = V + V (1) 


where V is the mean (time-averaged) velocity and v is 
the turbulent fluctuation component. The mean velocity, 
also known as the time-average, is given by: 

M 

(2) 


J= 1 


where M is the total number of samples at one boundary 
layer location. The root mean square value of the 
turbulent velocity fluctuation is: 


v 



E <Tj - V ) 2 (3) 

;=i 


and the local turbulence intensity is defined as: 


4*100 - 4 n 


— E (Yr Vf* loo 

mU 1 


( 4 ) 


For unsteady cases, the ensemble -averaged velocity, 
fluctuation velocity, and the turbulence intensity were 
calculated from the instantaneous velocity samples by: 


m) =- <m)> = IE w 


Nt 


/=i 


(5) 








( 6 ) 


Tu t (t)> = <Tu t {t)> = 


<m> 


xioo 


(7) 


where N= 100 is the total number of wake generator 
periods and M the number of samples taken per period. 
<Vi (tj)> is the reference ensemble averaged velocity for 
the particular boundary layer traverse. 

INTERMITTENCY ANALYSIS 

The intermittency distribution, which identifies 
whether the flow is laminar or turbulent inside the 
boundary layer, is calculated following the method of 
Hedley and Keffer [19]. Instantaneous velocities are 
used to identify this intermittency distribution. The 
instantaneous velocity is sensitized to increase its 
discriminatory capabilities between turbulent and non- 
turbulent parts of the signal. For this purpose, the 
multiplication of the first derivative of the velocity signal 
and the velocity signal is used for further analysis. This 
is called the detector function, S(t). Several other 
detector functions were used by Kowasznay et al. [21] 
and Antonia and Bradshaw [20]. 

( 8 ) 

Though sensitized detector function separates the 
turbulent and non-turbulent zones of the fluid, there is 
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still some overlap between two near the origin. The 
discrimination between the two zones of the flow will be 
ideal when the overlap between the two distributions is 
minimal or zero. To eliminate the disturbing effects of 
the velocity signal peaks, a smoothing procedure is 
applied to the S(t) signal: The mean value of ten 
consecutive S(t) values is calculated and the ten values 
are substituted by their mean value of S sm (t). 

After smoothing the detector function, a threshold 
level C is then applied to the smoothed detector function 
to distinguish between true turbulence and the signal 


noise. 


m 




when 

when 


SJf) 

SJf) 


C, 

c. 


( 9 ) 


After applying the threshold level to the detector 
function S(t), the result is a random square wave with 
0's representing the laminar case and 1 's representing the 
turbulent behavior of the boundary layer. A threshold 
level, C, of 1.2 is used for all the data on the suction 
surface. In the absence of length scales, this value is 
chosen from visual observations. Several other values of 
C are tested and little qualitative difference is seen in the 
intermittency di stribution during transition. Though the 
intermittency values vary with different values of C, the 
important parameters like start and end of transition are 
not effected by C. The resulting square wave after 
applying the threshold is ensemble- averaged to get the 
ensemble- averaged intermittency as follows: 


<YA)> = 

n J= 1 


( 10 ) 


Where n is the number of revolutions of the wake 
generator for which the data are collected. For time- 
averaged intermittency, <Y i (t i )> is integrated with 
respect to time to arrive at: 


— T 

Y=jf <Y ,(*,)><# 


( 11 ) 


t = 0 


Figure 3 show the processing of instantaneous 
velocities. 




Figure 3. Calculation of ensemble-averaged 
intermittency function from instantaneous 
velocities for Q=1.725 at y=0.720 mm. 


RESULTS AND DISCUSSION 

Time Averaged Velocity and Intermittency 
Distributions 

The effect of the wake frequency on the time- 
averaged velocity profiles and fluctuation velocity 
distribution are presented for one steady and two 
unsteady inlet flow conditions on the suction surface 
along 3 1 streamwise locations for the Reynolds number 
of 1 1 0,000. The steady state case serves as the reference 
configuration. After completing the velocity 
measurements, the boundary layer coordinates were 
transformed into a blade orthogonal coordinate system. 
Velocities at blade normal positions were obtained by 
interpolating their transformed values. The results 
showed almost no difference between the interpolated 
and non-interpolated velocity data. Experimental 
investigations were performed for three different values 
of Q = 0.0, 1.59, and 3.18. These values cover the 
reduced frequency range encountered in LPT-design 
and off-design operation conditions. 

The velocity fluctuation and intermittency 
distributions are shown in Figures 4 to 6 at 6 
representative streamwise locations for upstream region 
of the separation bubble, where the flow is attached. 
Upstream of the separation bubble at s/s c = 0.52 and also 
at its immediate proximity s/s o =0.588, the velocity 


7 


Copyright © 2005 by ASME 


distributions inside the boundary layer experience a 
slight decrease with increasing the reduced frequency. 
Inside the separation bubble at s/s 0 =0.705, a substantial 
influence of the wake frequency is observed. The higher 
wake frequency introduces a fluctuation kinetic energy 
into the boundary layer trying to reverse the separation 
tendency. As it can be seen from the velocity 
distribution profiles in Figure 4, the onset and the length 
of the separation bubble are not changed. However, 
there is a slight change of the bubble height. This shows 
that the wake flow does not have the capability to 
suppress the separation bubble. It only reduces the 
separation bubble height. In the downstream of the 
separation bubble, where the flow is fully reattached, 
s/s o =0.951, the impact of the wake on the boundary 
layer is reduced. This effect is clearly shown in the 
velocity distribution at s/s o =0.951. According to the 
previous investigations reported in [28] on a HP- 
turbine cascade, an increased wake frequency causes 
turbulence fluctuations to rise inside and outside the 
boundary layer. However, in the LPT- case with the 
boundary layer separation, once the boundary layer is 
re-attached and the velocity distribution assumes a fully 
turbulent profile, no major changes are observed neither 
in the velocity nor in the fluctuation distribution, Figure 
5. The time averaged intermittency distributions are 
shown in Figure 6 for several streamwise locations. 
Upstream of the separation bubble, the intermittency is 
very low through s/s 0 = 0.502. Close to the separation 
leading edge at s/s o =0.52 a sudden increase of 
intermittency is observed that indicates a change in flow 
state from a transitional to a fully turbulent state. As 
Figure 6 shows, inside the separation bubble at 
s/s o =0.705 intermittency changes drastically for all three 
frequency cases. An initial increase in intermittency is 
influenced by the shear layer responding with a slight 
decrease in intermittency. For a streamwise location of 
s/s o =0.705, this occurs at a lateral position ofy = 8 mm. 
Outside the bubble, the intermittency increases 
approaching a value close to unity. This indicatives the 
continuation of a transition process that started at the 
bubble leading edge above the shear layer that separates 
the separation zone with the external flow. A similar 
situation is observed further downstream at s/s o =0.767. 
Once the separation bubble is left behind ( s/s o =0.951), 
the effect of the shear layer is still felt. Moving from the 
wall, intermittency increases, reaches a peak at 
approximately y = 15 mm. It experiences a relaxation 
with a minimum value of <y(t)> = 0.5 which 
corresponds to the shear layer intermittency value. From 
this point on, the transition process along the shear layer 
determines the intermittency picture. It is worth noting 


that in all cases depicted in Figure 6, highest 
intermittency occurred in steady case with Q = 0. Inside 
the bubble, the wake effect seems not be strong, 
however, outside the bubble, the well known bee aiming 
effect makes itself felt in terms of reducing the 
intermittency. 

Ensemble- Averaged Intermittency Distribution 

The temporal-spatial contours of the ensemble- 
averaged intensity distribution at three different lateral 
positions for two reduced frequencies are presented in 
Figure 7. The wakes periodically disturb the boundary 
layer with the high turbulence intensity cores. As it is 
seen in Figure 7, the first three wakes are shown for 
better comparison of the effects of the impinging wake 
frequency. In these figures, the wakes with the highly 
vortical cores display intermittency values close to unity 
indicating the turbulent character of the boundary layer 
at the particular instant of time that the wake impinges 
on the surface. Intermittency is approximately equal to 
zero outside the wake region near the leading edge 
showing the non-turbulent behavior of the flow. The 
wakes represented by thin strips pass through the 
turbine blade channel and periodically switch the 
boundary layer from laminar to turbulent and vice versa 
depending on the presence of the wakes. At s/s o =0.52, 
the visibility of the wake is vanished due to the 
separation bubble. As explained earlier, the separation 
bubble starts at s/s o =0.52 and extends up to s/s o =0.746, 
thus occupying more than 24% of the suction surface 
and forming a massive separation. At s/s o =0.746, the 
intermittency field in Figure 7(a) shows the stagnant 
fluid region, which indicates the development of a 
transition and re-attachment. Increasing Q to 3.18 
causes an earlier mixing of the impinging wakes, which 
leads to a complete degeneration of the deterministic 
periodic flow into a stochastic turbulent flow. 

Relative Intermittency Distributions 

The intermittency distributions in Figure 7 clearly 
show the unsteady nature of the boundary layer 
transition. In this form, however, they cannot 
quantitatively describe the complex unsteady transition 
process. To establish the basic relations essential for a 
quantitative description of the unsteady boundary layer 
transition, we resort to the fundamental studies by 
Schobeiri and his co-workers ([14], [27], [30]) that deal 
with the physics of steady and unsteady wake 
development in a curved environment. These studies 
show that the turbulence structure of the steady and 
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Re=1 10,000, s/s o =0.502 


Re=1 10,000, s/s o =0.52 


Re=1 10,000, s/s o =0 .588 





Re=1 10,000, s/s o =0.705 



Re=1 10,000, s/s o =0.767 Re=1 10,000, s/s o =0.951 



Figure 4: Distribution of time-averaged velocity along the suction 
surface for steady case Q=0 (S R =°°) and unsteady cases Q=1.59 
(S r =160 mm) and Q=3.18 (S R =80 mm) at Re=l 10,000 


R©=1 10,000, s/s o =Q502 


Re=1 10,000, s/s o =0.52 


Re=1 10,000, s/s o =0 588 



Re=1 10,000, s/s o =0.705 Re=1 10,000, s/s o =0.767 Re=1 10,000, s/s o =0 .951 





Figure 5: Distribution of time-averaged fluctuation rms velocity along 
the suction surface for steady case Q=0 (S R =°°) and unsteady cases 
Q=1.59 (S r =160 mm) and Q=3.18 (S R =80 mm) at Re=l 10,000 
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Re=1 10,000, s/s 0 =0.951 



Figure 6: Distribution of time -averaged intermittency along the suction 
surface for steady case Q=0 (S R =°°) and unsteady cases Q=1.59 (S R =160 
mm) and Q=3.18 (S R =80 mm) at Re=l 10,000 


unsteady wake flow is determined by the wake defect, 
which is a Gaussian function. Following the above 
studies, we define a dimensionless parameter: 


< 12 > 


that relates the passing time, t, of a wake impinging on 
the plate surface with the wake passing velocity in the 
lateral direction, U w , and the intermittency width, b. 
The latter is directly related to the wake width 
introduced by Schobeiri and his co-workers [30]. In an 
analogous way to find the defect function, we define the 
relative intermittency, T, as: 


r = in 

<Y - < Y,(0 > n m, 


In the above equation, <y^)> is the time dependent 

ensemble-averaged intermittency function, which 
determines the transitional nature of an unsteady 
boundary layer. The maximum intermittency 


<V ft"! > exhibits the time dependent ensemble- 

‘ i'V |/max 

averaged intermittency value inside the wake vortical 
core. Finally, the minimum intermittency 


< yA ) 1 


. > represents 

mm 


the ensemble- averaged 


intermittency values outside the wake vortical core. 
Figure 8 exhibits the maximum and minimum ensemble- 
averaged intermittency inside and outside the wake 
vortical core. 

A representative relative intermittency function, T, 
is shown in Figure 9 (a, b, c, d) for a frequency value of 
Q=1.59 at a lateral distances from the blade surface of 
y =0.858, 0.996, 5.3, and 9.3 mm, with the 
dimensionless longitudinal distance s/s 0 as a parameter. 
The above lateral distances are representative for 
intermittency distributions inside, within and outside the 
separation bubble over the entire suction surface. The 
symbols represent the experimental data. As seen, for 
the reduced frequency of Q=1 .59, the measured relative 
intermittency functions follow very closely a Gaussian 
distribution, given by: 


r = 


(14) 
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Q=3.18, y= 1.341 mm 







<*t)> 

1.00 

0.95 

0.88 

0.82 

0.76 


0.63 
0.57 
— 0.50 
0.44 
0.38 
0.32 
0.25 
—I 0.19 
— | 0.13 
0.06 
0.00 


Figure 7. Ensemble-averaged intermittency factor in the temporal-spatial domain at different y positions 
for Q=1.59 ( S r =160 mm), and Q=3.18 ( S R =80 mm) 
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and represented by the solid curve. Here, £ is the non- 
dimensionalized lateral length scale from Eq. (12). 
Using this function as a generally valid intermittency 
relationship for unsteady wake flows, the intermittency 
function <y i (t i )> is completely determined if additional 
information about the minimum and maximum 
intermittency functions <Y / (/j) max > and < y i (^ y ) rn - ri > 

are available. Figure 11 shows the time-averaged 
intermittency distribution for one steady or no rod case 
and two unsteady cases on the suction surface, as a 
function of s/s 0 at different normal positions from the 
blade. Upstream of the leading edge of the separation 
bubble, the time averaged intermittency is determined by 
the laminar nature of the boundary layer, which exerts 
a strong damping effect on the impinging wake 
fluctuations as extensively discussed in [29]. 
Approaching the bubble leading edge a steep increase in 
intermittency indicate a strong turbulent fluctuation 
within the separation bubble. This exactly corresponds 
to the findings plotted in Figure 5. Close to the wall ( y 
~ 0. 1,0. 720) the intermittency peak is embedded in the 
separation bubble as shown in Figure 11. Moving 
toward shear layer causes an increases the value of the 
peak intermittency. At s/s 0 = 0.705, where the separation 
bubble height reaches its maximum, the intermittency 
approaches its minimum and increases again to reach the 
second maximum. The change of the intermittency state 
is reflected in Figure 10, which was extensively 
discussed in [29]. 



Figure 8. Maximum and minimum intermittencies at 
y=0.1 mm and s/s o =0.383 


The distributions of <y i (t l ) ms > ar, d <Y ? (^ / ) rn ^> 

in streamwise direction are plotted in Figure 12 for Q = 
1.59 and 3.18 on the suction surface. The distribution 
of< Yi(^) max > corres P on ^ s to th e condition when the 


wake with its high turbulence intensity core impinges on 
the plate surface. Once the wake has passed over the 
surface, the same streamwise location is exposed to a 
low turbulence intensity flow regime with an 
intermittency state of < y i(Pmin> > where no wake 


present. As shown in Figure 12, outside the separation 
zone the minimum intermittency < y tenc ^ s to 


follow the course of steady (no wake) intermittency 
distribution. The final state of <y (t\ . > does not 

* iV i/ min 


approach the fully turbulent value of 1.0 due to the 
wake calming effect. 


CONCLUSIONS 


A detailed experimental and theoretical study ofthe 
behavior ofthe separation bubble on the suction surface 
of a highly loaded LPT blade under periodic unsteady 
wake flow is presented. The measurements were carried 
utilizing a custom designed hot-wire probe. One steady 
and two different unsteady inlet wake flow conditions 
with the corresponding passing frequencies, the wake 
velocity and the turbulence intensities were investigated 
by utilizing a large-scale, subsonic research facility. Two 
dimensional wakes were generated by cylindrical rods 
attached to two parallel timing belts performing a 
continuous translational motion in front of the turbine 
cascade. The results of the unsteady boundary layer 
measurements and the intermittency analysis were 
presented in the ensemble -averaged, and contour plot 
forms. The analysis of the boundary layer experimental 
data with the flow separation, confirms the universal 
character of the relative intermittency function which is 
described by a Gausssian function. 

The minimum intermittency factor, <Y min >, 
represented the intermittency, when the boundary layer 
is subjected to the wake external region. While upstream 
of the separation bubble, the minimum intermittency 
< Y min > , resembles the distributions found in [1], [2] and 
[3], within the separation bubble a steep drop followed 
by a moderate increase dictates the intermittency 
picture. 
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(C) 


£1=1 .59, y=5.3 mm 


(d) 


£1=1 .59, y=9.3 mm 



Figure 9. Relative intermittency as a function of s/s 0 for unsteady frequency of 
Q=1.59 (S r =160 mm) at (a) y=0.858 mm, (b) y=0.996 mm, (c) y=5.3 mm, and (d) 
y=9.3 mm at Re=l 10,000 


(a) £1=1.59, tft=0.25 (b) £1=1 .59, t/x=0.50 



Figure 10. Ensemble -averaged velocity contours along the suction surface for different s/s 0 with time t/T 
as parameter for Q=1.59 (S R =160 mm), Re=l 10,000 (time-averaged separation bubble for Q=1.59 
marked red) 
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(e) y=2.1 mm (f) y=1 0.1 mm 



Figure 11. Time-averaged intermittency as a function of s/s c at different lateral positions for steady case 
Q=0 (S R =°°) and unsteady cases Q=1.59 (S R =160 mm) and Q=3.18 (S R =80 mm) at Re=l 10,000 
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(a) Q=1.59, y=0.765 mm (b) Q=3.18, y=0.765 mm 



Figure 12. Maximum, minimum and time-averaged intermittency as a function of s/s G at different lateral 
positions for steady case Q=0 (S R =°°) and unsteady cases Q=1.59 (S R =160 mm) and Q=3.18 (S R =80 
mm) at Re=l 10,000 
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A more detailed picture of the intermittency 
behavior inside the laminar, separation, and the turbulent 
region is given by the time averaged intermittency. 
Upstream of the separation bubble the intermittency is 
determined by the laminar nature of the boundary layer, 
which exerts a strong damping effect on the impinging 
wake fluctuations as extensively discussed in [29]. 
Approaching the bubble leading edge a steep increase in 
intermittency indicates a strong turbulent fluctuation 
within the separation bubble. Close to the wall, the 
intermittency peak is embedded in the separation bubble. 
Moving toward the shear layer causes an increase in the 
value of the peak intermittency. The intermittency 
approaches a minimum, where the separation bubble 
height reaches its maximum. 

UNCERTAINTY ANALYSIS 

The Kline and McClintock [31] uncertainty analysis 
method was used to determine the uncertainty in the 
velocity after calibration and data reduction for the 
single- wire probe. In addition, the uncertainty in the 
heat transfer measurements was also determined. The 
Kline and McClintock method determines the 
uncertainty with a 95% confidence level. The 
uncertainty in the velocity for the single-wire probe after 
the data reduction is given in Table 3. As shown, the 
uncertainty in the velocity increases as the flow velocity 
decreases. This is due to the pneumatic pressure 
transducer having a large uncertainty during calibration. 
Table2: Uncertainty in velocity measurement for hot- 
wire probe. 


U (m/s) 

3 

5 

12 

coU/U (%) rod 

5.78 

2.41 

1.40 
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